APS/123-QED 



Indirect exchange coupling between localized magnetic moments in carbon nanotubes: 

a dynamic approach 

A. T. CostE0 and R. B. Muniz 
Instituto de Fisica, Universidade Federal Fluminense, 24210-346 Niteroi, RJ, Brazil 

M. S. Ferreirsd] 

School of Physics, Trinity College Dublin, Dublin 2, Ireland 
(Dated: February 29, 2008) 

Magnetic moments dilutely dispersed in a metallic host tend to be coupled through the conduction 
electrons of the metal. This indirect exchange coupling, known to occur for a variety of magnetic 
materials embedded in several different metallic structures, is of rather long range, especially for 
low-dimensional structures like carbon nanotubes. Motivated by recent claims that the indirect 
coupling between magnetic moments in precessional motion has a much longer range than its static 
counterpart, here we consider how magnetic atoms adsorbed to the walls of a metallic nanotube 
respond to a time-dependent perturbation that induces their magnetic moments to precess. By 
calculating the frequency-dependent spin susceptibility we are able to identify resonant peaks whose 
respective widths provide information about the dynamic aspect of the indirect exchange coupling. 
We show that by departing from a purely static representation to another in which the moments are 
allowed to precess, we change from what is already considered a long range interaction to another 
whose range is far superior. In other words, localized magnetic moments embedded in a metallic 
structure can feel each other's presence more easily when they are set in precessional motion. We 
argue that such an effect can have useful applications leading to large-scale spintronics devices. 

PACS numbers: 



I. INTRODUCTION 

The demand for faster computers has transformed the 
field of magnetoelectronics from a purely academic field 
of research into one of the most promising scientific ar- 
eas regarding potential technological applications. When 
combined with systems of reduced dimensionality like 
thin films, nanowires and molecular structures, this field 
of research is expected to pave the way to the produc- 
tion of non-volatile computer memories, efficient mag- 
netic sensors and magnetic materials with enhanced in- 
formation storage capacity, to name but a few. 

Among low-dimensional structures, carbon nanotubes 
are potentially useful for magneto-transport applications. 
The ability to produce sizeable changes in the conduc- 
tance of a nanotube due to an applied magnetic field is 
one of the driving forces in the research of magnetic prop- 
erties of carbon-based structures [TJ [2]. In this case, nan- 
otubes must interact with magnetic foreign objects that 
lift the intrinsic spin balance of a non-magnetic material 
like carbon. Substrates [3j H] , substitutional impurities [5 1, 
adsorbed atoms j6[ [7] and nanoparticles[8] are some of 
the different magnetic foreign objects that can interact 
with carbon nanotubes. Among those, transition-metal 
magnetic adatoms have been reported to produce notice- 
able changes in the spin-dependent electronic structure 
of carbon nanotubes [51 [B]. Furthermore, defect-induced 
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magnetic moments in carbon-based materials have been 
reported as another way to manipulate the magnetic re- 
sponse of these low-dimensional systems |S] • 

Motivated by the interaction of carbon nanotubes with 
foreign magnetic objects, we have recently investigated 
the problem of magnetic coupling between localized mag- 
netic moments embedded in metallic carbon nanotubes 
[TOl [TTJ [12] . The direct exchange interaction between 
magnetic moments is known to decay abruptly as a func- 
tion of their separation since this interaction requires a 
finite overlap between the wave functions that surround 
the respective magnetic objects. When the moments are 
not too close together there is no direct overlap between 
their wave functions and the only way a magnetic cou- 
pling can arise is if it is mediated by the conduction elec- 
trons of the metallic host. This so-called indirect ex- 
change coupling (IEC) between localized magnetic mo- 
ments mediated by the conduction electrons of metallic 
hosts often plays a central role in determining the overall 
magnetic order in dilute magnetically-doped metals. 

Put another way, the IEC measures how far apart lo- 
calized magnetic moments embedded in a metal can feel 
each other's presence. Metallic nanotubes fall within the 
category of materials that may display long ranged cou- 
plings and have been predicted to display some very in- 
teresting features when doped with magnetic obiects fTUl 
[TTJdZj. More specifically, in Reference [TO] we have con- 
sidered how two magnetic adatoms attached to the walls 
of a nanotube and separated by a distance D are mutu- 
ally coupled. Surprisingly, we found that the coupling 
lacks the typical oscillations observed in indirect inter- 
action of magnetic objects embedded in 3-dimensional 
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metallic hosts. We have explained this lack of oscilla- 
tions as a commensurability effect responsible for hiding 
the truly oscillatory character of the IEC. Moreover, we 
have shown that the IEC in nanotubes decays very slowly 
giving this interaction a very long range. Reference 
deals with the competition between direct and indirect 
exchange couplings that occurs when moments are not 
too far apart and points to the possible existence of non- 
collinear magnetic order in dimerized impurities. Finally, 
reference [12] shows that the long range character of the 
IEC depends on the exact location of the magnetic mo- 
ments relatively to the underlying atomic structure of the 
carbon nanotubc. 



All the features mentioned above are associated with 
the magnetic coupling in a purely static configuration, 
that is, without accounting for any dynamic aspect of the 
spin degrees of freedom. Recently, there has been a grad- 
ual shift of interest from the static nature of the magnetic 
interaction towards dynamic aspects of the magnetism of 
low-dimensional systems, primarily driven by the tech- 
nological importance of fast magnetization switchings in 
memory devices. Regarding the magnetic coupling, re- 
cent reports indicate the existence of an indirect interac- 
tion between precessing magnetic moments separated by 
a non-magnetic metallic multi-layered material |13[ 114] . 
In particular, precessing magnetic moments that are cou- 
pled despite being very far apart suggests that the range 
of the so-called dynamic coupling always exceeds that of 
its static counterpart by a considerable amount. Bear- 
ing in mind that the static IEC in nanotubes is already 
very long ranged [THl HU H2] , it is only natural that we 
test whether the magnetic coupling can indeed reach fur- 
ther when magnetic moments embedded in nanotubes 
are allowed to precess. This is the goal of the current 
manuscript. For the sake of distinction to earlier descrip- 
tions of the dynamic coupling, it is worth mentioning 
that here we adopt a full quantum mechanical formalism 
that includes contributions not present in standard semi- 
classical formalisms, like the one originally carried out in 
Reference [15] . 



The sequence adopted in this manuscript is the follow- 
ing. Next section presents the geometry of the system 
we will be working with, together with a definition of the 
corresponding Hamiltonian that describes its electronic 
structure. The spin susceptibility is a key quantity in 
this work and is also introduced in the following section. 
After a physical interpretation is given for the resonant 
peaks of the diagonal matrix elements of the susceptibil- 
ity, we shift our interest towards its off-diagonal matrix 
elements which will be then used to quantify the dynamic 
coupling throughout this paper. Section |rTT] shows the re- 
sults of the dynamic coupling between magnetic adatoms 
for different cases, followed by a brief section with con- 
clusions. 




FIG. 1: (a) Schematic diagram of two magnetic adatoms, la- 
beled A and B, adsorbed to the walls of a carbon nanotube. 
They are separated by a distance D. (b) Possible positions 
for the magnetic atoms relatively to the underlying hexagonal 
lattice. On the left, the magnetic atom is adsorbed on top of 
a single carbon atom, also referred to as the atop position. 
In the diagram, the adatom has been slightly shifted for clar- 
ity. The middle diagram shows the bridge position in which 
the adatom is above the line joining two nearest-neighbour 
carbons. Finally, on the right, in the central position, the 
adatom is above the centre of the underlying hexagon. 



II. MODEL DESCRIPTION 

We consider two magnetic atoms, labelled A and B, 
adsorbed to the walls of an infinitely long metallic car- 
bon nanotube (armchair) and schematically represented 
in Fig. [IJl The precise location of these so-called 
adatoms is one of the three possibilities shown in Fig. 
[TJd. The first possibility, depicted by the left diagram, 
corresponds to the adatom being immediately above one 
carbon atom[T5]. Alternatively, it can lie symmetrically 
above the bonding line between two neighbouring car- 
bon atoms (middle) or it can be right above the centre of 
a hexagon composed of six carbon atoms (right). They 
are labeled as atop, bridge and central positions, respec- 
tively. It has been recently shown that the range of the 
static IEC depends on the precise location of the adatom 
relatively to the underlying hexagonal '.at 1 ice 12 . It is 
therefore instructive to investigate if and how the dy- 
namic coupling depends on the adatoms location. 

We assume that the magnetic moments of the indi- 
vidual adatoms A and B arc initially parallel. We refer 
to this initial setup as the ferromagnetic configuration, 
which can be the result of a spontaneous parallel align- 
ment [10] , depending on the adatoms separation, or can 
be induced by an externally applied uniform magnetic 
field H that drives all the moments to be mutually paral- 
lel. Let us now assume that an additional time-dependent 
magnetic field localized on atom A is applied perpendic- 
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ularly to the original magnetization direction inducing 
the moment of adatom A to precess with a frequency u>. 
Since we are interested in the indirect coupling between 
the moments A and B, we must investigate the effect 
that this induced precession on the former will have on 
the latter. In other words, by calculating the induced 
motion of the magnetization on adatom B driven by the 
time-dependent field applied on adatom A, we have a di- 
rect way of assessing how strong the indirect coupling 
between the moments is. This is given by the spin sus- 
ceptibility x( w ); a matrix whose elements reflect how the 
spin degrees of freedom of the system respond to a mag- 
netic excitation. 

To calculate the spin susceptibility one needs the 
Hamiltonian operator describing the electronic structure 
of the system, which is defined by H = Hmt + Ha + 

H B +Vc + H z , where H NT = E^v^T^ £ j> is the 
Hamiltonian of the individual nanotube. The opera- 
tors c' m a and c m ,o- are the creation and annihilation op- 
erators, respectively, for a 7r-electron of spin a = f , J, 
located at site m of the nanotube. The sum over j 
and j' are over nearest neighbours only and the pa- 
rameter 7 stands for the electronic hopping 7 = 2.7eV. 
Ha and Hb are Hubbard-likc Hamiltonians defined as 
Hp = Eo-^^.ff^" + ^ ™A,crn A ^), where (3 = A,B 
and h x = c x c x is the number operator. The parameters 
tp and Up represent the on-site potential and the intra- 
atomic electronic interaction on the adatom j3, respec- 
tively. The operator V c = T"£) f)<7 (c] )0 .CA,<T + c\ a c Lc! ) + 

T * J2e> crC^l' <A B >° + a^i' ,<?) describes the interaction 
between the adatoms and the nanotube, where r plays 
the role of an electronic hopping between them. The in- 
dices £ and £' label the nanotube sites that are coupled 
to adatoms A and B, respectively. Finally, Hz plays the 
role of a Zeeman term in the Hamiltonian and is propor- 
tional to the applied magnetic field. 



ih r + °° 

3&,» = ^ / du'f{J){[Gl m {J)~G}l n ^')]G, 



In the equation above n (ui) is the Fourier trans- 
form of the time-dependent single-particle Green func- 
tion Gm,n(t) f° r an electron with spin a defined as 

^,„W = -^0W({w(i),4, CT }> ■ (3) 

In what follows, we show that one way of describ- 
ing how precessing magnetic moments are aware of each 
other's presence is by looking at the frequency-dependent 
spin susceptibility. Different matrix elements of this 
quantity reflect distinct features of spin excitations, some 
of which will be subsequently discussed. 



We choose to describe the electronic structure of the 
system by the single-band tight-binding model, which is 
justified by the fact that the main features of the IEC 
are predominantly determined by the extended electrons 
of the host, in this case the nanotube. As a matter 
of fact, the electronic structure of carbon nanotubes is 
known to be well reproduced by such a model. With 
transition-metal atoms in mind, the adatoms are also de- 
scribed by a single orbital representing a 5-fold degener- 
ate (i-band with the appropriate occupation to represent 
typical magnetic materials. The parameter r is selected 
based on comparisons with ab-initio evaluations for the 
static IEC [TT] and the parameter Up was assumed to 
be 1.677 throughout this paper. Despite the simplified 
model for the electronic structure, our approach is by 
no means restricted to such a simple case. The results 
here obtained can be easily extended to a multi-orbital 
description but bring no qualitative difference. 

The time-depedent transverse susceptibility is defined 
as X m,n(t) = ~i®(t)([S+(t),S-(0)}), where 0(aj) is the 
heaviside step function, and S + and S~ are the spin rais- 
ing and lowering operators, respectively. The indices m 
and n refer to the locations where the field is applied and 
where the response is measured, respectively. In our case, 
inducing the moment of adatom B to precess and mea- 
suring the corresponding fluctuations of the magnetiza- 
tion on adatom A is described by XA,B{t). After Fourier 
transforming to frequency domain, the susceptibility can 
be further manipulated and rewritten in matrix form as 

xM = [i + L>x°M]- 1 x°M , (l) 

where U is a diagonal matrix in which U a and Ub are the 
sole elements and \° i s the Hartree-Fock susceptibility 
whose matrix elements are 



>' + w) + [G^ n {J) - G% m (u/)] G^JiJ + u)} . (2) 
I 

III. RESULTS 

A few characteristic plots of the frequency-dependent 
spin susceptibility are displayed in Fig(2] whose features 
require further analysis. The figure shows the imagi- 
nary part of matrix elements of the spin susceptibility 
x(cj). These quantities are directly associated with how 
strongly the system responds to the time-dependent mag- 
netic field. In other words, large values of x(w) indicate 
a greater sensitivity of the system to excitations of fre- 
quency to. In particular, peaks in the w-dependent sus- 
ceptibility reflect the existence of resonant frequencies 
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FIG. 2: Plots of -Imxm,n(w) (in arbitrary units) as a func- 
tion of hui (in units of 7), where m and n label the respective 
matrix element of the susceptibility. All diagrams are for a 
(5,5) armchair nanotube and for magnetic atoms described by 
a d-band occupation of 6 electrons (e.g. Fe). (a) m — n = A 
corresponds to the diagonal element of the susceptibility cal- 
culated for a single magnetic atom (A) adsorbed on the cen- 
tral position, (b) Similar matrix element (m = n = A) for 
two magnetic atoms (labeled A and B) separated by D = 3a, 
where (c) Off-diagonal matrix element of the susceptibility 
for the same system as in (b). The difference is in the matrix 
element, now defined by m — A and n = B. 



and whose respective widths characterize the lifetime of 
the spin excitation. In the example depicted on the top 
panel of Fig. [2j the imaginary part of xa,a{w) corre- 
sponding to an armchair nanotube with only a single 
magnetic adatom (labeled A) has a clearly distinctive 
peak. On the middle panel, Ithxa,a{^>) is calculated for 
two identical adatoms (A and B) a distance D = 3a 
apart along the longitudinal axis of the nanotube, where 
a is the lattice parameter of the graphcnc lattice. Two 
separate peaks are clearly identifiable and correspond to 
resonant frequencies associated with the different normal 
modes which the system can oscillate with, namely the 
acoustic and the optical modes. In the former the mo- 
ments of both atoms A and B oscillate in phase but are 
out of phase in the latter mode. 

Bearing in mind that the peak width is inversely pro- 
portional to the characteristic lifetime of a pulsed exci- 
tation, it is possible to determine from panels (a) and 
(b) of Fig. [2] whether or not the magnetic moments 
are dynamically coupled. In the case of a single adatom 
the peak is visibly broader than the corresponding peak 
(acoustic mode) for the case of two adatoms. This in- 
dicates that the precession of atom A is shorter-lived 
when the adatom is isolated. It is worth highlighting 
that since there is no spin-orbit interaction in the Hamil- 
tonian, there are no dissipative terms capable of lessening 
the total angular momentum of the system. This means 



that the finiteness of the precessional motion lifetime has 
nothing to do with dissipation but must be related to 
the way in which the magnetic moments interact with 
the embedding medium. 

As a matter of fact, we interpret this difference in 
the excitation lifetimes by suggesting that the precession 
of the magnetic moment interacts with the surrounding 
conduction electrons of the nanotube host and induces 
an angular momentum current, often referred to as spin 
current |13l 114] , that propagates away from the process- 
ing magnetic moment. In the case of a single magnetic 
atom, this spin current travels away from the atom with 
no probability of ever returning to the location where it 
was originally generated. When two magnetic atoms are 
present, the spin current emitted by one atom is capable 
of exciting the precession motion of the other and vice- 
versa. In this case, the spin current emanating from the 
atoms is not entirely lost and can be repeatedly scattered 
back and forth, which explains why the acoustic mode on 
Fig|2j) is narrower than that of Figj2k. 

Variations in the width of the acoustic resonant peak 
are therefore indicative of a dynamic coupling between 
coexisting magnetic atoms and can in principle be used to 
probe the existence of this indirect interaction. Fig(3^i il- 
lustrates how this can be done by plotting Af2/f2o for two 
magnetic adatoms separated by a distance D , where Af2 
is the width of the resonant peak in the acoustic mode 
and fi is the precise value of the resonant frequency. 
The width Af2 is defined as the energy range around 
Q,q within which the susceptibility exceeds \a,a(^o) /2. 
A horizontal dashed line highlights the corresponding di- 
mcnsionless quantity for a single impurity. Despite a sud- 
den discontinuity at D — 4a, reminiscent of the standard 
oscillatory behaviour found in this type of interaction [T4"] . 
the general trend is an increase in the width of the acous- 
tic peak towards the limiting value displayed by a single 
isolated magnetic atom. 

Besides providing a simple physical picture that ex- 
plains the observed results in the susceptibility calcula- 
tions, the main reason why the dynamic coupling is an- 
alyzed in terms of the widths of the resonating peaks is 
because this is traditionally how ferromagnetic resonance 
experiments probe excitations of this nature. Neverthe- 
less, this is unfortunately not the most direct way for 
evaluating the dynamic coupling between magnetic im- 
purities and an alternative definition is desirable. 

Rather than inferring about the coupling by variations 
in the width of the resonant peaks, the imaginary part 
of the off-diagonal matrix element of the spin suscepti- 
bility connecting the two magnetic atoms A and B is 
also a representative quantity for this type of interac- 
tion. 1xcixa,b(w) describes how the moment at atom 
A responds to an induced excitation of frequency uj at 
atom B and should therefore reflect the behaviour seen 
in Fig(3^i in a more direct way. In fact, in spite of not 
having a clear proportionality relation with the curve in 
Figj3^i, Fig. [3Jd displays — JmXA,B (^0) as a function of D 
and shows a similar trend to the one reported above in 
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FIG. 3: Two different ways of inferring about the dynamic 
coupling. On (a) the width AQ of the resonant peak (centred 
at the acoustic-mode frequency Qq) of —lm\A,A is plotted as 
a function of the adatom separation D (in units of a). When 
compared with the same quantity obtained for a nanotube 
with only a single magnetic adatom (represented by the hor- 
izontal dashed line), we see that spin excitations of adatoms 
that are relatively close together live longer than those for 
which the adatoms are isolated. On (b), similar trends are 
shown for —lm\A,B (in arbitrary units), which appears as 
an alternative quantity displaying the same features. Once 
again, in both cases the adatoms are located on the central 
positions and are chosen with the same d-band occupation as 
those of Fig. [2] 



terms of the resonating peak widths, that is, of a coupling 
that decreases with the atomic separation D. Notice that 
even the discontinuity at D — 4a is reproduced. It is 
thus evident that by associating the dynamic coupling 
with the imaginary part of xa,b{^o) we obtain a much 
simpler way of assessing how this interaction depends on 
D. This simplicity is invaluable given the aforementioned 
motivation of testing whether the range of the dynamic 



FIG. 4: — lmxA,B(Qo) (in arbitrary units) is plotted as a 
function of the adatom separation D, this time for adatoms on 
the atop position. fio corresponds to the resonant frequency 
of the acoustic mode. Solid line displays the calculated values 
for the dynamic coupling. Starting from D = 10a, the dashed 
line indicates how the coupling should fall if it was to follow 
the same rate of decay as the static IEC. 



interaction is indeed superior to that of its static coun- 
terpart. 

Before discussing how the dynamic coupling that ap- 
pears in Fig. [3]3 scales with the adatom separation, it is 
worth mentioning that the static IEC has been reported 
to decay unusually fast (with approximately l/D 5 ) when 
magnetic adatoms are located on the central positions of 
metallic nanotubesjT^]. This fast decaying rate for the 
static IEC was shown to result from interference effects 
characteristic of the underlying hexagonal lattice of nan- 
otubes. The dynamic coupling represented by the imagi- 
nary part of xa.b is depicted by the solid line of Fig. [SJa 
and decays considerably more slowly than l/D 5 . As a 
matter of fact, starting from D = a, the dashed line rep- 
resents how the dynamic coupling would decay if it was 
to fall with the same rate as that of the static coupling. 
Regarding the solid-line curve of Fig. despite a rapid 
decay (of above 2 orders of magnitude) between D = a 
and D = 10a, this is still far slower than what one would 
expect from the purely static interaction. 

For adatoms on both the atop and bridge locations, the 
static IEC has been shown to decay very slowly with the 
adatom separation D. In fact, the observed l/D decaying 
rate for these two cases is what characterizes the static 
IEC as a long-range interaction [15] . When the dynamic 
aspects of the indirect coupling are included, the range of 
this interaction is extended even further, as can be seen in 
Figs. [4] and [5] Dashed lines on both graphs indicate the 
l/D behaviour expected for the static coupling whereas 
the solid lines show the dynamic coupling through the 
imaginary part of the off-diagonal matrix element of the 
spin susceptibility. In both cases, a pronounced oscilla- 
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FIG. 5: — ImxA,s(fio) (i n arbitrary units) is plotted as a 
function of the adatom separation D, this time for adatoms on 
the bridge position. f2o corresponds to the resonant frequency 
of the acoustic mode. Solid line displays the calculated values 
for the dynamic coupling. Starting from D — 9a, the dashed 
line indicates how the coupling should fall if it was to follow 
the same rate of decay as the static IEC. 

tory behaviour is evident but most remarkably, the range 
of the dynamic coupling is far superior to that of the 
static IEC. In summary, what was already considered a 
long range interaction appears to reach even longer dis- 
tances when the magnetic moments on the walls of a 
metallic nanotube are put in precessional motion. 

IV. CONCLUSIONS 

The direct conclusion from the results presented above 
is that no matter how long ranged the static IEC is, 



the dynamic version of this indirect interaction seems 
to be of superior range. This is in agreement with 
earlier semi-classically-based predictions that a dynamic 
coupling may exist even in the absence of its static 
counterpart |13l 114] . Here, based on a quantum mechan- 
ical approach describing the electronic structure of the 
system, we are able to confirm this earlier prediction and, 
moreover, conclude that magnetic moments attached to 
the walls of a metallic nanotube can feel each other's 
presence further apart when they are in motion. The fact 
that setting the moments in motion expands the range 
of the indirect interaction is quite remarkable, bearing in 
mind that the static IEC is already very long ranged due 
to its reduced dimensionality. 

Being this dynamic interaction of such long range, it 
is likely that this could be exploited as a communica- 
tion channel between distant parts of spintronic devices. 
Schemes for synchronizing large sets of oscillators have 
also been devised based on the dynamic coupling between 
magnetic moments |16l [T7]. In particular, the ability to 
excite in-phase oscillations to a magnetic moment that is 
far apart from the region where the inductive magnetic 
field is originally produced may enable control of very 
large arrays of coupled spin-transfer devices opening the 
road to the construction of macroscopically-sized devices 
with the above discussed quantum features. 



Acknowledgments 

The authors would like to thank D. F. Kirwan for his 
assistance with some of the figures in this manuscript. M. 
S. F. acknowledges the financial support of Science Foun- 
dation Ireland. A. T. C acknowledges support received 
from the brazilian research council, CNPq. 



[1] K. Tsukagoshi, B. W. Alphenaar, and H. Ago, Nature 

401, 572 (1999). 
[2] L. E. Hueso, J. M. Pruneda, V. Ferrari, G. Burnel, J. 

P. Valdes-Herrera, B. D. Simons, P. B. Littlewood, E. 

Artacho, A. Fert, and N. D. Mathur, Nature 445, 410 

(2007) 

[3] M. S. Ferreira and S. Sanvito, Phys. Rev. B 69, 035407 
(2004). 

[4] O. Cespedes, M. S. Ferreira, S. Sanvito, M. Kociak, and 

J. M. D. Coey, J. Phys.: Condens. Matter 16, 155 (2004). 
[5] S. B. Fagan, R. Mota, R. J. Baierle, A. J. R. da Silva, 

and A. Fazzio, Mater. Charact. 50, 183 (2003). 
[6] S. B. Fagan, R. Mota, A. J. R. da Silva, and A. Fazzio, 

Phys. Rev. B 67, 205414(2003). 
[7] S. B. Fagan, R. Mota, R. J. Baierle, A. J. R. da Silva, 

and A. Fazzio, Physica B 340, 982 (2003). 
[8] C. K. Yang, J. Zhao, and J. P. Lu, Phys. Rev. Lett. 90, 



257203 (2003) 

[9] Yuchen Ma, A. S. Foster, A. V. Krasheninnikov, and R. 
M. Nieminen, Phys. Rev. B 72, 205416 (2005) 

[10] A. T. Costa, D. F. Kirwan, and M. S. Ferreira, Phys. 
Rev. B 72, 085402 (2005) 

[11] A. T. Costa, C. G. Rocha and M. S. Ferreira, Phys. Rev. 
B 76, 085401 (2007) 

[12] D. F. Kirwan, C. G. Rocha, A. T. Costa and M. S. Fer- 
reira, Phys. Rev. B 77, 085432 (2008) 

[13] B. Heinrich, Y. Tserkovnyak, G. Woltersdorf, A. Brataas, 
R. Urban and G. E. W. Bauer, Phys. Rev. Lett. 90, 
187601 (2003) 

[14] A. T. Costa, R. B. Muniz and D. L. Mills, Phys. Rev. B 

73, 054426 (2006) 
[15] In the diagram, the adatom has been slightly shifted for 

clarity. 

[16] S. Kaka, M. Pufall, W. Rippard, T. Silva, S. Russek, and 



7 



J. Katine, Nature 437, 7057 (2005) 
[17] J. Grollicr, V. Cros, and A. Fert, Phys. Rev. B 76, 060409 
(2006) 



